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ABSTRACT 



Aims. We present parallaxes of 1 1 mid-to-late T dwarfs observed in the UKIRT Infrared Deep Sky Survey. We use these results to 
test the reliability of model predictions in magnitude-color space, determine a magnitude-spectral type calibration, and, estimate a 
bolometric luminosity and effective temperature range for the targets. 

Methods. We used observations from the UKIRT WFCAM instrument pipeline processed at the Cambridge Astronomical Survey 
Unit. The parallaxes and proper motions of the sample were calculated using standard procedures. The bolometric luminosity was 
estimated using near- and mid-infrared observations with two different methods. The corresponding effective temperature ranges were 
found adopting a large age-radius range. 

Results. We show the models are unable to predict the colors of the latest T dwarfs indicating the incompleteness of model opacities 
for NH3, CH4 and H2 as the temperature declines. We report the effective temperature ranges obtained. 

Key words. Astrometry - Stars: low-mass, brown dwarfs, fundamental parameters, distances 



1. Introduction 



Since the first discovery of a T dwarf by Nakaiima et al .I (Il995h 
understanding the atmospheric processes and the role of chemi- 
cal composition for such low temperature objects has been very 
challenging. This remains an important goal, particularly given 
that we consider these objects to be the link between stars and gi- 
ant exoplanets, therefore, their properties offer insights into for- 
mation and evolution of planetary systems. 

One of the fundamental parameters required to determine 
the physical properties and to constrain theoretical models of 
celestial objects is distance. Until now, only a few late T 
dwarfs have known distances: W olf940 B (lHarrington & Dahn 



98 0]), HD3651B dPerrvman et al"1ll997l), 2MASS J0415-093 5 
( Vrba et alJl2004h . 2MASS .109 39-2448 dBurgasser et alT 2008). 
ULAS J003402.77-005206.7 (ISmart et alJ 120101) . Thus the 
theoretical models have been constrained b y the earlier T 
dwarf s. The UKIRT Infrared Deep Sky Survey (lLawrence et alJ 
I2007L hereafte r UKIDSS) is revealing large numbers of 



new T dwarfs riBurningham et all 1201 (A fLodieu et al.1 l2009allb|: 



Burningham et al. 2008; Pinfield et al. 2008; Chiu et al. 2008 



Lodieu et alJ 120071 ; Kendall et al. 2007) and in particular cool 



T8/T9 dwarfs. An important follow up for the coolest T dwarfs 
observed in the UK IDSS is the deter mination of their par- 
allax. As shown in Smart et a D (f20Tol hereafter SJL10) for 
ULAS J003402.77-005206.7 using the UKIDSS discovery im- 
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age it is possible to reduce the time required to determine a pre- 
liminary parallax and here we present 10 additional objects. 

In this contribution all targets will be referred to by the 
discovery acronym and right ascension short format hence 
ULAS J003402.77-005206.7 becomes ULAS 0034, the full- 
names are given in Table Q] In Section 2 we describe the ob- 
servations and procedures and in Section 3 we report the results 
for the 1 1 targets. In Section 4 we compare these results to cur- 
rent models and in Section 5 we calculate L^ ; and T e ff of the 
objects in the sample. Finally in Section 6 we discuss the results 
obtained. 



2. Observations and Reduction Procedures 

The observations for the parallax determination began in 2007 
and the target list of 11 objects was drawn from confirmed 
T dwarfs observed at that time in the UKIDSS Large Area 
Survey. In Table Q] we list th e objects along with MKO mag- 
nitudes dTokunaga et alj2 002) and near-infrared (hereafter NIR) 
spectral types. The N IR spectral type classi fication follows the 
scheme described in Burga sser et alJ (120061) . The only object 
with an already published parallax was SDSS 0207 which was 
include d because the US NO parallax had a large relative error 
(~28%. IVrbaetai]l2004l) . 

The procedures for observing, image treatment and parallax 
determination follow those described in SJL10. In that paper the 
parallax observations of ULAS 0034 maintained the same target 
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position as the discovery image and were aligned using a simple 
linear transformation. Assuming the astrometric distortion pat- 
tern does not change, the use of the discovery image as part of 
the parallax solution allows a much shorter dedicated observa- 
tional campaign to obtain a precise parallax. For CFBDS 0059, 
SDSS 0207, ULAS 0948 and ULAS 2239 this was not possi- 
ble because the target was very close to a chip edge in the first 
UKIDSS image. The WFCAM astrometric distortion is signif- 
icant, so moving the reference frame on the focal plane results 
in poorer astrometric transformations (i.e. larger residuals in the 
solutions) and so in this situation those frames that are signifi- 
cantly offset are given lower weight. 

Initially all observations were made with the same total ex- 
posure time of 400 seconds. For this exposure, in the case of the 
anonymous stars in the field of ULAS 0034, the centroiding pre- 
cision is a constant 20mas until around /=18.4 and deteriorates 
quickly to 60mas at 7=19.2 (Figure 2 in SJL10). To ensure that 
centroiding precision is optimal even in poor observing condi- 
tions we have increased the exposure times to 710s for the faint 
targets (ULAS 0034, 0948, 1 150, 1315, 2239 and CFBDS 0059). 
The precision of the final solutions for the faint targets reflects 
the poor quality of some of the earlier 400 second observations. 

In figure Q] we plot the solution of ULAS 0901 (top panel) 
and 1315 (bottom panel) which are objects of similar distances 
and observational history but with magnitudes of 17.90 and 
18.86 respectively, e.g. straddling the precision borderline, to 
show high and low quality solutions. 

3. Results 

The astrometric parameters derived for the 1 1 targets are in Table 
[2] For each one we report: target name, position (J2000), number 
of reference stars, number of observations used, absolute par- 
allax, proper motion components, tangential velocity, the time 
span covered by the observations and the relative-to-absolute 
parallax correction applied. The proper motion of the targets have 
all been brought to an absolute system using the galaxies in the 
field. 

The relative errors for most of the targets are less than 10%. 
Observations are continuing and we hope to be able to reduce 
them to 5% by the end of the campaign. Significant exceptions 
are ULAS 1150, 1315 and 2239 which are also amongst the 
faintest objects under study, hence degraded the most by the 
borderline signal-to-noise in the early observations. In addition, 
parallax observations of ULAS 2239 were shifted from the dis- 
covery image position lowering the weight of the first point. 

In Table [3] we compare the astrometric distances obtained 
here with the estimated ones given in the discovery paper of each 
object. The estimations were made with different techniques and 
the reader is referred to the discovery papers (see footnote to 
TableQ]) for further details. The discovery distance range is usu- 
ally much larger than our measured one and is overestimated for 
the late-Ts. This comparison underlines the need for measured 
parallaxes that are model independent. 

4. Model Fitting 

In Figs. |2] and [3] we present four color - absolute magnitude 
diagrams for a sample of L and T dwarfs with My > 12 and 
> 12 respectively. The 11 targets are plotted as filled circles. 
Parallaxes and magnitudes of the other objects are from two 
on-line archives: 

1) The L/T dwarf archive maintained by S. K. Leggett 
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Fig. 1. Observations for the targets ULAS 0901 (top panel) and 
1315 (bottom panel). The highest point is the discovery image. 
The observational history and distance of these two targets are 
similar but they differ by a magnitude in apparent brightness. 
The solution shows the effect of low signal-to-noise observations 
in the beginning of the parallax sequence. 



(http://staff.gemini.edu/~sleggett/2OlO41ho1_tab.txt, hereafter 
Leggett archive): This archive contains a compendium of 225 
objects with MKO YJHKL'M' magnitudes and IRAC [3.55], 
[4.49], [5.73] and [7.87] magnitudes. Where used these objects 
are plotted as filled squares. 

2) The M/L/T dwarf archive (www.dwarfarchives.org, hereafter 
Dwarf archive): This is an on line compendium of all published 
L/T dwarfs and selected M dwarfs. As of 10/01/2010 there were 
752 L & T dwarfs, reporting JHK magnitudes, parallaxes and 
proper motions. Where used these objects are plotted as filled 
triangles. 

To make Figs. [2] and [3] clearer we omitted the error bars on 
the literature objects. An indication of the typical uncertainty on 
these points is given by the black cross above the legend. 
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Table 1. Infrared magnitudes and spectral types of the 1 1 targets. 



Full Name 


ZAB 


Y 


J 


H 


K 


Sp. Type 


Refs (D,P,T) 


ULAS J003402.77-005206.7 


22.11 ± 


0.05 


18 


90 = 


= 


10 


18 


15 ± 





03 


18 


49 


+ 





04 


18.48 


±0.05 


T9 


1,1,6 


CFBDS J005910.90-01 1401.3 


21.93 + 


0.05 


18 


82 h 


- 


02 


18 


06 ± 





03 


18 


27 


+ 





05 


18.63 


±0.05 


T9 


2,2,6 


SDSS J020742.48+000056.2 


20.11 + 


0.60 


17 


94 h 


= 


03 


16 


75 + 





01 


16 


79 


+ 





04 


16.71 


±0.05 


T4.5 


3,(4,7),8 


ULAS J082707. 67-020408. 2 






18 


29 i 


= 


05 


17 


19 ± 





02 


17 


44 


+ 





05 


17.52 


±0.11 


T5.5 


4,4,4 


ULAS J0901 16.23-030635.0 






18 


82 h 


= 


05 


17 


90 ± 





04 


18 


46 


+ 





13 


> 1 


8.21 


T7.5 


4,4,4 


ULAS J094806.06+064805.0 






20 


03 i 


= 


14 


18 


85 + 





07 


19 


46 


+ 





22 


> 18.62 


T7 


4,4,4 


ULAS J101821.78+072547.1 






18 


90 i 


= 


08 


17 


71 + 





04 


17 


87 


+ 





07 


18.12 


±0.17 


T5 


4,4,4 


ULAS Jl 15038.79+094942.8 


22.44 ± 


0.10 


19 


92 h 


= 


08 


18 


68 ± 





04 


19 


23 


+ 





06 


19.06 


±0.05 


T6.5p 


5,5,5 


ULAS J131508.42+082627.4 


22.82 ± 


0.10 


20 


00 4 


= 


08 


18 


86 ± 





04 


19 


50 


+ 





10 


19.60 


±0.12 


T7.5 


5,5,5 


ULAS J133553.45+1 13005.2 


22.04 + 


0.10 


18 


81 i 


= 


04 


17 


90 ± 





01 


18 


25 


+ 





01 


18.28 


±0.03 


T9 


6,6,6 


ULAS J223955. 76+003252.6 






19 


94: 


= 


17 


18 


85 ± 





05 


19 


10 


+ 





10 


18.88 


±0.06 


T5.5 


4,4,4 



YJHK are in the MKO Vega photometric system, while z is in AB syste m. The uncertainty i n th e spectral type i s ±0.5 . 

References (D = D isc overy, P = Photome try T = Spectral Type): 1- [Warren et all d2007h 2- IDelorme et all d2008h 3- iGeballe et ail d2002h 4- 
iLodieu et alj d2007h 5- lPinfield et"ai] d2008h 6-lBurningha m et alj d2008l) 7- lKnapp et alJd2004h 8- lBurgasser et alJd2006ft 



Table 2. Parallaxes and proper motions of the 1 1 targets. 



Target 


a 


S 


N*, obs 


rf bs ± <r n 


abs , 
H-a — u fa 


abs , 
H-6 — "M 


v te „ 


Time span 


COR 




(h:m:s) 


(d:m:s) 




(mas) 


(mas/y) 


(mas/y) 


(km/s) 


(years) 


(mas) 


ULAS 0034 


0:34:02.7 


- 0:52:07.8 


135, 17 


78.0 ± 3.6 


-18.5 ±3.2 


-363.3 ± 3.6 


21.7 ± 1.0 


3.81 


1.24 


CFBDS 0059 


0:59:10.9 


- 1:14:01.4 


70, 13 


108.2 ±5.0 


878.8 ± 8.4 


50.5 ± 4.8 


38.6 ± 1.8 


2.69 


1.19 


SDSS 0207 


2:07:42.9 


+ 0:00:56.0 


47, 17 


29.3 ± 4.0 


158.8 ±3.0 


-14.3 ±3.9 


25.8 + 3.6 


3.70 


1.37 


ULAS 0827 


8:27:07.6 


- 2:04:08.4 


418, 17 


26.0 ±3.1 


26.8 ± 2.7 


-108.9 ±2.3 


20.5 ± 2.5 


3.94 


0.87 


ULAS 0901 


9:01:16.2 


- 3:06:35.4 


241, 19 


62.6 ± 2.6 


-38.6 ± 2.3 


-261.2 ± 2.8 


20.0 ± 0.8 


3.90 


1.00 


ULAS 0948 


9:48:06.1 


+ 6:48:04.5 


152, 15 


27.2 ± 4.2 


199.4 ± 7.0 


-273.9 ± 6.2 


59.1 ±9.3 


1.58 


0.98 


ULAS 1018 


10:18:21.7 


+ 7:25:46.8 


198, 14 


25.0 ± 2.0 


-183.7 ±2.6 


-15.1 ±3.1 


34.9 ± 2.8 


2.00 


1.04 


ULAS 1150 


11:50:38.7 


+ 9:49:42.8 


105, 10 


16.8 ±7.5 


-107.6 ± 17.1 


-31.9 ±4.5 


31.7 + 14.9 


2.91 


1.09 


ULAS 1315 


13:15:08.4 


+ 8:26:27.0 


213, 11 


42.8 ± 7.7 


-60.2 ± 8.3 


-95.8 ± 10.0 


12.5 ± 2.5 


3.07 


0.95 


ULAS 1335 


13:35:53.4 


+ 11:30:05.1 


196, 8 


96.7 ± 3.2 


-196.9 ±4.9 


-201.0 ± 6.3 


13.8 ±0.5 


2.18 


0.95 


ULAS 2239 


22:39:55.7 


+ 0:32:52.7 


120, 15 


10.4 ± 5.2 


125.3 ± 5.4 


-108.4 ± 5.2 


75.7 ± 37.8 


2.95 


0.96 



In the fourth column we report the number of reference stars (N*) and the number of observations (obs). In the ninth column we report the time 
span covered by the observations and in the last one the relative-to-absolute parallax correction (COR). 



Table 3. The 1-sigma distance range obtained here compared to 
those estimated in the discovery papers. 



Name 


Astrometric 


Discovery 


Discovery 




distance (pc) 


distance (pc) 


reference 


ULAS 0034 


12.2-13.4 


14-24 


1 


CFBDS 0059 


8.8-9.6 


8-18 


2 


ULAS 0827 


33.8-43.0 


24-39 


3 


ULAS 0901 


15.3-16.7 


21-33 


3 


ULAS 0948 


31.1-42.5 


38-60 


3 


ULAS 1018 


36.8-43.2 


33-52 


3 


ULAS 1150 


32.9-86.1 


42-60 


4 


ULAS 1315 


19.2-27.6 


34-48 


4 


ULAS 1335 


10.0-10.6 


8-12 


5 


ULAS 2239 


48.0-144.2 


52-83 


3 



References: 1 - IWarren et al.l d2007h 2- IDelorme et alj d2008|) 3- 
iLodieu et alj d2007h 4- iPinfield et alj d2008h 5- iBurningham eTai] 
d2008h 



All magnitudes are in the MKO system and preferentially 
taken from the Leggett archive as they are measured directly 
in that system. The majority of the infrared magnitudes in the 
Dwarf archive are in the 2MASS system and when needed we 
use the relations in Stephens & Leggettj (|2004|) to convert to the 
MKO system. 

Colored lines i n Figs . [2] and [3] are the model predictions 
by iBurrows et aD (120061 hereafter BSH06) and lAUardetalJ 



(I2003L 120071 120091 hereafter BTSettl09). BSH06 tracks covers 
the temperature range 700-2000 K, with log[g]=4.5,5.0,5.5 and 
[Fe/H]=-0.5,0,+0.5. BTSettl09 covers the range 500-780 K, with 
log[g]=4.5,5.0,5.25 and [Fe/H]=-0.2,0,+0.2. For a given metal- 
licity (indicated in the plots by a given color) log[g] increases 
from left to right in J-H, while it increases from right to left in 
J-K. 

Model colors and magnitudes were obtained by convolv- 
ing the theoretical spectra with the UKIDSS filter profiles 
dHewett et al.l2 006) to calculate fluxes. We interpolate the model 
spectra with a spline to have the same binning as the filter pro- 
files, apply the selected profile, and integrate to obtain the to- 
tal flux. The integrated fluxes were then converted into an abso- 
lute magnitude using as a zero point a Vega spectrum treated the 
same way. 

The absolute magnitudes plotted in Figs. [2] and [3] require 
from the models the flux at 10 pc. The BSH06 models sup- 
ply the flux at the surface of the object and at 10 pc, the 
latter calculated ass uming the radius-log(g)-T e /y relation from 
IBurrows et ail dl997l) . The BTSettl09 models provide only the 
flux at the surface o f the object, to find t he flux at lOpc we as- 
sume the radii from Baraffe et al.l d2003l) corresponding to the 
model log(g) and effective temperature. 

In Figs. [2] and [3] we note that the BTSettl09 tracks for high 
and low metallicity are swapped in the different color spaces. 
In J-H space tracks for high metallicity are bluer than the low 
metallicity ones, while in J-K they are redder. This follows the 
trend suggested by BSH06 tracks for higher temperature, and 
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can be seen also in the models of lSaumon & Marlevl (l2008). Due 
to this swap, the predictions for the T9 dwarfs are incompatible 
in the two color spaces, i.e. they predict high gravity - low metal- 
licity in J-H and low gravity - high metallicity in J-K. Moreover, 
the T9s are redder than predicted by the theoretical models, espe- 
cially in J-K. This failure may be due to an incorrect prediction 
of the flux emitted in the H and K-band resulting from known 
opacity deficiencies in modeling the collision-induced absorp- 
tion of H2 and the wing of the K I doublet resonance, and the 
lack of an appropriate lis t of absorption line s of CH4 and NH3 at 
such low temperatures (Leggett et al. 2010l). 

The predictions for earlier objects are consistent in the two 
color spaces, but in Fig. 2 they tend to slightly underestimate the 
magnitudes of the T6.5-T7s. Further discussion on the models 
predictions and individual objects is in Section 6. 

In Fig. H we present three absolute magnitude - IR spec- 
tral type diagrams. The IR spectral type cl assification follows 
the scheme described in iGeballe et al. (2002) for L dwarfs and 
the scheme described in BurgasseretaL (2006) for T dwarfs. 
The o ver plotted curves are polynomial fits derived in iLiu et ail 
(2006), with the dotted line representing the polynomial ob- 
tained excluding from the fit all the known and possible bina- 
ries, while the dashed line is the one obtained excluding only 
the known binaries. Possible binaries were selected by Liu et al. 
based on their relatively high T e f f compared to objects of simi- 
lar spectral type in the iGolimowski et alj d2004t) measurements. 
One of the possible binaries indicated b yLiuetal. (SPSS J1 021- 
0304) was confirmed as a binary by Burga sser et al.l d2006l) . so 
now is plotted as a known binary. Following the convention used 
before, filled circles are the 11 targets and squares are objects 
from the Leggett archive, while blue objects are known binaries 
and green objects are possible binaries. The trend for T8.5 and 
T9 is an extrapolation of the Liu et al polynomial, since their 
sample consisted of objects between L0 and T8. 

The 1 1 targets indicate a steeper trend in the sequence be- 
yond T8 than the extrapolation of the Liu et al. polynomial. In 
Fig. |4] the red lines are our 4-th order polynomials fit to the data 
including the new objects presented here, but excluding the 14 
objects without magnitu de naturally in the MKO system. This 
choice was made as the Ste phens & Leggett! d2004l) transforma- 
tions from the 2MASS to the MKO system were derived using 
hotter objects and employing them may introduce systematic er- 
rors in the resulting fit of the cooler objects. In Table [4] we list 
the coefficients and errors of the fit. The difference between the 
two fits with and without possible binaries is reduced compared 
to Liu et al., probably because of the smaller statistical weight 
of the possible binaries in this study (5 over a total sample of 61 
objects) with respect to Liu et al. (6 over 43). As the sample of 
late T dwarfs is still small, identification of possible binaries in 
this region may be incomplete. 

We also tested polynomials from 3rd to 10th order but after 
4th order there was no significant improvement in the sigma of 
the fit. 

5. Luminosity and Effective Temperature. 

Next we calculate an effective temperature range for our targets. 
To do this we used the classical Stefan-Boltzmann law: 

L bol = ^(rR 2 T A eff (1) 

hence to calculate the temperature we need to know the radius 
and the bolometric luminosity of each object. 

To determine a radius, we consider the models of 
Baraffeet~ai1 d2003l) we see that the radius of T dwarfs decreases 



rapidly when the object is very young but after 0.5 Gyr it is very 
constant. Therefore if we can constrain the age of our objects 
we can use these models to constrain the radius. In Fig. [5] we 
report the galactic velocity components (U, V,W) of the 1 1 tar- 
gets. The components were calculated from the proper motions 
in Table [2] assuming a radial velocity range of +80/-80 km/s. 
Given this large range in V ra d< the uncertainty in the proper mo- 
tion becomes negligible, so we ignore it in this calculation. The 
over plotted box is the locus o f very young objects (age < 0.1 
Gyr,|Zuckerman & Song 2004J) w hile the ellipso id is the locus of 
young disk objects (age < 0.5 Gyr, Eg genii 9691) . Even given this 
very conservative range in possible radial velocities no objects 
fall in the boxed area or in the ellipsoid area in all three compo- 
nents. We therefore conclude that all obje cts are older tha n 0.5 
Gyr. For the 0.5-10 Gyr models from iBaraffe et all d2003l) we 
find the radius has a range of 1.2 - 0.8 Rj up (in a temperature 
range of 500 - 2000 K). This i s also consistent with the radii pre- 
dicted by Bur rows et alj d 1997b models for objects of the same 
age and temperature. 

To calculate the bolometric flux (and hence the luminosity) 
we combined the available measured spectr um of each object 
dBurningham et al]l2008HLeggett et al.ll2009h flux calibrated us- 
ing UKIDSS YJHK photometry, with the model spectra. For 
ULAS 0034 and 1335 we have the spectrum in the near- and 
mid-infrared (hereafter MIR) region, while for the other objects 
we only have the NIR portion. To find the bolometric flux we use 
model spectra to estimate the flux where we do not have obser- 
vations: at short wavelengths {A < 1.0 /urn), and, for ULAS 0034 
and 1335, in the region between the near- and mid-infrared spec- 
trum (2.4 //m < A < 7.5 fim) or, for the other targets, in the entire 
portion from 2.4 yum to 15 //m. The flux emitted beyond 15 /im 
was estimated assuming a Rayleigh-Jeans tail. 

The models used are the alr eady mentioned BTSettl0 9 
and the non-equilibrium models by Huben v~& Burrowsl (|2007), 
covering the temperature range from 700 to 1900 K, for 
log[g]=4.5,5.0 and 5.5, assuming values of K, z =10 2 ,10 4 and 
10 6 cm 2 s _1 (eddy diffusion coefficient), fo r different speed of 
the CO/CH4 reaction (for further details see lHubenv & Burrowsl 
l2007t ISaumon et al]|2006l and reference therein). 

For each object we took the models in a wide range of 
temperature (+200 K around the temperature predict ed by the 
temperature-type relation given bv lStephens eUd"] |2()0971). grav- 
ity (log[g] from 4.5 to 5.5) and metallicity ([Fe/H] from -0.2 to 
+0.2) and we scaled them using the available magnitudes (listed 
in Tables Q] and |5j- For each model we took the average scal- 
ing factor obtained and we compared it with the range given by 
the known distance and the radius range adopted (i.e. the square 
of 0.8 Ry„ ; ,/distance + 3cr - the square of 1.2 Ry 1(p /distance + 
3cr). We discarded the models whose average scaling factor was 
out of this pseudo-3cr range. We joined each one of the remain- 
ing spectra with the measured one and we calculate the result- 
ing bolometric flux, the luminosity and hence the temperature 
range (corresponding to the radius range). Then we compared 
the temperature obtained with the one associated with the model 
employed. We discarded those models whose temperature dif- 
fered by more than 100 K from the temperature range obtained. 
Finally we assumed the mean flux given by the remaining mod- 
els as our final estimation and hence we calculate the luminosity 
and the temperature range. 

The uncertainty in the flux is given by the spread in values 
plus a 3% uncertainty in the magnitudes used to calibrate the 



1 Except for the T9s, since we don't have theoretical spectra for tem- 
peratures lower than 500 K. 
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Fig. 2. Color-magnitude diagrams for a sample of L and T dwarfs. The colored lines are theoretical tracks from BSH06 and 
BTSettl09, for different gravities and metallicities. For each metallicity, the gravity increases from left to right, assuming the values 
4.5,5.0,5.5 for BSH06 tracks and 4.5,5.0,5.25 for BTSettl09. The 11 targets presented here are plotted as filled circles with their 
associated error bars, the Dwarf archive objects are plotted as triangles and the Leggett archive objects as squares. All magnitudes 
in the MKO system. The cross above the legend indicates the typical uncertainty on the literature objects. 



observed spectra and to scale the model ones. The uncertainty in 
the luminosity and the temperature is the result of the standard 
propagation of the errors on the flux and the distance, ignoring 
the uncertainty in the radius, given the wide range adopted. 



The results are shown in Table [6] In the first column we indi- 
cate the target short name, in the second one its spectral type, 
in the third the temperature esti mated using the tem perature- 
spectral type relation given by Stephens et al. (2009), in the 
fourth the range of models employed, in the fifth the range of 
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Fig. 3. Same as Fig. [2]but with the color J-K instead of J-H. For each metallicity, the gravity increases from right to left, assuming 
the values 4.5,5.0,5.5 for BSH06 tracks and 4.5,5.0,5.25 for BTSettl09. 



models kept after the first step (so after comparing the scaling 
factors), in the sixth the range of models kept after the second 
step (so after comparing the temperatures), in the seventh our 
assumed bolometric flux, in the eighth the associated luminosity 
and in the last the temperature range obtained. 

We note that the use of MIR magnitudes and spectra in- 
creases our ability to constrain the object temperature. For 
ULAS 0034 and 1335 the bolometric flux is well constrained 



and the width of the temperature range obtained (150-200 K) is 
mainly due to the radius. Using the NIR spectrum and photom- 
etry only, the uncertainty on the flux increases and the tempera- 
ture range obtained doubles. 



To have an idea of the eventual systematics, we also tested 
the technique described in Cushing et al. (2008), e.g. to fit the 
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IR Spectral type 



Fig. 4. /, H and K-band absolute magnitude as a function of IR spectral type. Labeled objects are those with new parallaxes presented 
here. Blue squared points are know n to be unresolve d binaries, while green squar ed points are possibl e unresolved binaries. The over 
plotted lines are polynomial fit by iLiu et all |2006) based on data tabulated by Knap p et al.l d2004l) . The dotted line was obtained 
excluding all the known and possible binaries, the dashed one excluding only the known binaries. The red lines are our new 
polynomial fits. 
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Table 4. Coefficients of the polynomial fit. 



Magnitude 


ao 




ai a 2 




a 3 






34 




Excluding known and possible binaries 


Mj 


11.140 ± 


0.052 


0.556 ±0.032 3.74 ± 0.62 xlO" 2 


-9.49 ± 


0.45 xlO 


-s 


3.69 


±0.11 xlO 


-4 


M H 


10.486 ± 


0.051 


0.590 ± 0.031 -3.80 ± 0.61 xl0~ 3 


-4.35 ± 


0.44 xl0~ 


-3 


2.15 


±0.11 xlO- 


4 


M K 


10.180 ± 


0.056 


0.375 ± 0.035 2.31 ± 0.67 x 10" 2 


-5.03 ± 


0.49 xlO" 


-3 


2.11 


±0.12 xlO- 


4 


Excluding known binaries 


Mj 


10.957 ± 


0.053 


0.736 ± 0.032 -5.53 ± 0.61 xl0~ J 


-6.31 ± 


0.44 xlO" 


-J 


2.96 ±0.11 xlO" 


-4 


M H 


10.292 ± 


0.052 


0.778 ±0.031 -4.77 ± 0.60 xl0~ 2 


-1.13 ± 


0.43 xlO" 


-3 


1.41 


±0.10x10 


-4 


M K 


10.005 ± 


0.056 


0.548 ±0.034 -1.84 ± 0.64 xlO" 2 


-1.93 ± 


0.46 xlO- 


-3 


1.38 


±0.11 xlO- 


4 



The fits are defined as in lLiuetal.lJ2006h : M xx = Z%a a, x (SpT)' where XX indicates J, H or AT-band magnitude and the spectral types are defined 
following the convention SpT=l for LI, SpT=9 for L9, SpT=10 for TO etc. The fit is valid for spectral types from LO to T9. Infrared spectral types 
are used for both L and T dwarfs. 

Table 6. Fluxes, luminosities and temperatures of the sample obtained scaling the model spectra using the measured magnitudes, 
e.g. the first method discussed in Section 5. 



object 


Sp. 


est. 


models 


1st step 


2nd step 


F^o/ 




L/L G 


T e // 


range (cr) 




Type 


Teff (K) 


used (K) 


(K) 


(K) 


(erg s _I cnT 2 


) 






(K) 


ULAS 0034 


T9 


500 


500:700 


540:700 


540:660 


2.34±0.05xl0- 


-13 


1.13±0.06xl0- 6 


535 


- 660 (35) 


CFBDS 0059 


T9 


500 


500:700 


500:620 


500:620 


2.85±0.48xl0- 


-13 


7.20±1.25xl0- 7 


480 


- 590 (55) 


SDSS 0207 


T4.5 


1130 


900:1400 


1000:1400 


1100:1200 


4.98±0.26xl0- 


-13 


1.72±0.25xl0- 5 


1060- 


1300(110) 


ULAS 0827 


T5.5 


1070 


900:1300 


1000:1300 


1000:1100 


2.87±0.09xl0- 


-13 


1.21±0.15xl0- 5 


970- 


1190 (95) 


ULAS 0901 


T7.5 


830 


600:1000 


600:700 


600:700 


1.88±0.15xl0- 


-14 


1.42±0.13xl0- 6 


570 


- 700 (50) 


ULAS 0948 


T7 


910 


700:1100 


700:800 


700:800 


6.52±0.44xl0- 


-14 


2.61±0.44xl0~ 6 


660 


- 810 (80) 


ULAS 1018 


T5 


1100 


900:1300 


900:1100 


900:1000 


1.80±0.43xl0 


-13 


8.54±0.70xl0- 6 


890- 


1090 (75) 


ULAS 1150 


T6.5p 


980 


800:1200 


800:1200 


900:1000 


6.76±0.40xl0 


-14 


7.12±3.22xl0~ 6 


850- 


1040 (160) 


ULAS 1315 


T7.5 


830 


600:1000 


600:740 


600:660 


7.47±0.72xl0- 


-14 


1.21±0.25xl0- 6 


545 


- 670 (70) 


ULAS 1335 


T9 


500 


500:700 


540:660 


540:660 


3.42±0.09xl0- 


-13 


1.09±0.06xl0- 6 


530 


- 650 (35) 


ULAS 2239 


T5.5 


1070 


900:1300 


900:1300 


1100:1200 


6.00±0.09xl0- 


-14 


1.65±0.82xl0- 5 


1050- 


1280 (200) 



Table 5. IRAC m agnitudes of ULAS 0034, SPSS 0207 
and ULAS 1335 JWarren et alj 120071: iPatten et al.l l2006t 
Burningham et al . 2008) used to scale the model spectra. 





ULAS 0034 


SDSS 0207 


ULAS 1335 


3.55 //m 
4.49 /mi 
5.73 fim 
7.87 /mi 


16.28 ± 0.49 
14.49 ± 0.43 
14.82 ± 0.44 
13.91 ± 0.42 


15.59 ± 0.06 
14.98 ± 0.05 
14.67 ± 0.20 
14.17 ± 0.19 


15.96 ±0.48 
13.91 ± 0.42 
14.34 ± 0.43 
13.37 ±0.41 



object spectrum using the model spectra. The best fit spectrum 
was selected as the one that minimize: 

where n is the number of data pixels, f, is the measured flux in 
the i-th spectral interval, Q is the scaling factor (R/d) 2 , F^, is 
the k-th model flux and cr, is the error in the measured spectrum. 
The weight associated to each bin (w,) is the extension of the 
bin itself (A/l) as suggested by Cushing et al. The scaling factor 
can be provided by the fit, however, since we know the distance 
to the dwarf we consider fixed values of Q, assuming again the 
radius range previously indicated. We selected the best fit models 
for the two extreme configurations, i.e. 0.5 Gyr-1.2 Rj up and 10 
Gyr-0.8 R Jup . 

Since we don't have an associated noise spectrum for 
ULAS 0827, 0948, 1018, 1150, 1315 and CFBDS 0059 we min- 
imize: 

G * = L Wi — ?j— (3) 

;=1 V yji J 



The uncertainty in the extremes takes into account half of 
the models grid spacing based on fitting (for further details see 
Cushing et al.), plus an additional percentage on the flux due to 
the incomplete spectral coverage. This percentage was estimated 
comparing the measured IRAC magnitudes of ULAS 0034 and 
1335 with the model's predicted ones. The differences between 
measured and model magnitudes gives an average uncertainty 
of ~70% on the calculated flux between 2.5 and 7.5 //m. For 
ULAS 0034 and 1335 in this interval there is -40% of the to- 
tal emergent flux, so we obtain a relative error 0> ;,. c /=.7x.4=.28. 
This implies an additional uncertainty of ~7% on the temper- 
ature. For the other 9 objects, the uncertainty in the tempera- 
ture was estimated extending the relative sigma calculated for 
ULAS 0034 and 1335 to the uncovered part of the flux (-60%, 
that results in an uncertainty cr Fre i-.lx.6=A2). 

The choice of the weight function is arbit rary. Different 
choice s ca n lead to different resu lts, as seen by ICushing et al.l 
(12008b and lStephens etal.1 (120091) . Given this, we prefer the re- 
sults obtained with the first method described. The values ob- 
tained with this spectral technique are summarized in Table [7] 
They are largely consistent with the ones obtained scaling the 
model spectra. 

We also performed a completely model independent flux cal- 
culation for ULAS 0034 and 1335, that have a completely mea- 
sured spectral energy distribution. We determined the bolometric 
flux emitted integrating the measured spectrum between 1 and 
2.5 /mi, then adding the flux emitted between 2.5 and 7.5 fim 
calculated using the IRAC magnitudes (assuming a constant flux 
distribution over the passbancQ), finally we integrated the mea- 

2 Any error that this assumption may introduce would be negligible, 
since the IRAC passbands are tight. 



F. Marocco et al.: Properties of 1 1 T dwarfs 



9 



S o 




Table 7. Temperatures obtained fitting the observed spectrum, 
e.g. the second method discussed in Section 5. 



Fig. 5. The galactic velocity components U, V and W obtained 
from the proper motions in Table [2] assuming a V ra j range of 
+80/-80 km/s. The triangle indicates the +80 km/s extreme, the 
diamond indicates the -80 km/s extreme. The overplotted black 
box is the locus of young stars dZuckerman & Song1l2004 age< 
0.1 G yr), the ellipsoid is the locus of young disk stars (Eggen 
1969, age< 0.5 Gyr). 



Name 


Sp. Type 


T-// 


range (cr) 








(K) 


ULAS 0034 


T9 


500 


- 580 (40) 


CFBDS 0059 


T9 


500 


- 540 (60) 


SDSS 0207 


T4.5 


1000 - 


1200 (120) 


ULAS 0827 


T5.5 


1100- 


1200 (130) 


ULAS 0901 


T7.5 


620 


- 720 (75) 


ULAS 0948 


T7 


700 


- 800 (90) 


ULAS 1018 


T5 


900- 


1100(120) 


ULAS 1150 


T6.5p 


800- 


1000(110) 


ULAS 1315 


T7.5 


540 


- 620 (65) 


ULAS 1335 


T9 


500 


- 600 (40) 


ULAS 2239 


T5.5 


1100- 


1300(135) 



IR spectral type 

Fig. 6. Temperature ranges plotted as a function of the spectral 
type. Uncertainties in each extreme point of the range are in 
Table [6] The uncertainty in spectral type is half subtype. The 
over plotted dotted l ine is the effec t ive te mperature-infrared type 
relation derived by Ste phens et al.l d2009t) . 



sured 7.5 - 15 fim spectrum. The flux emitted beyond 15 fim was 
estimated assuming a Rayleigh-Jeans tail. The results obtained 
for ULAS 0034 and 1335 with this approach are consistent with 
the one obtained with the other techniques described within the 
uncertainty quoted in Table [6] 



6. Discussion 

In Fig. [6] we present a T e yy-spectral type diagram of the 1 1 T- 
dwarfs of the sample, for each of which we plot the temperature 
range displayed in Table[6] When needed, objects have been off- 
set by +0.1 in spectral type, to avoid overlaps. Over plotted for 
comparison we h ave the effective tempe rature - infrared type re- 
lation derived bv lStephens et al.l (|2009). All the ranges are con- 
sistent with the relation except for ULAS 0948, 0901 and 1315, 
that are cooler than predicted. 

We now discuss individual objects with temperature values 
or indications of peculiarity in the literature. 

6.1. ULAS 0034 

There are several estim ations of the t emper ature of this object. 
In the discovery paper, IWarren et al.l (|2007) estimate a conser- 
vative range of 600 < T e /y < 700 K using the measured NIR 
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spectrum and the near- and mid-infrared photometry with a grid 
of solar-abundance BTSettl model, calibrated using the param- 
eters of 2MASS 0415, and a linear fit of T e// vs. H-[4.5] color 
for hotter stars; iDelorme et al 

and the same grid described in Warren et al. I (l2007h determine 
T eff * 670 K; in lLeggett etail d2009l) the pubhshed range is 550 
< T e jf < 600, obtained fitting the measured NIR and MIR spec- 
trum with model spectra bv lSaumon & Marlevi (l2()08): finally, in 
SJL 10 using the parallax re ported in this paper and the model fit- 
ting in Leggett et al.l d2009l) . using a model radius of 0.1 1 R© (as 
implied by the spectral fits constraints on gravity), derive a bolo- 
metric luminosity of L/L = 1.10+0. 01 xl0~ 6 , consistent with the 
temperature range 550-600 K. In this work, using NIR and MIR 
spectrum, the MKO NIR photometry and the IRAC MIR pho- 
tometry, we find a luminosity o f L/L o =1.13±0.06x l0~ 6 , consis- 
tent with the value obtained by Leggett et al.l d2009l) . The differ- 
ence between the temperature range derived here and by Leggett 
et al. is due to the different technique used: the Leggett et al. 
spectral fitting result is similar to that derived spectrally here, 
given in Table [7] 



6.2. CFBDS 0059 



trast with this finding) but we must poi nt out the known d egen- 
eracy between gravity and metallicity dKnapp et al.l l2004). This 
degeneracy could be the reason for the discrepancy between our 
2008) using t he NIR spectrum result and Lodieu et al. 



In the discovery paper dDelorme et al.ll2008h of this object they 
estimate T e yy~620 K comparing the spectral indi ces with a so- 
lar me tallicity grid of BTSettl model spectra. In Legge tt et al.l 
(2009() the technique adopted is the same described in Sec. 
6.1, and the range obtained is 550<T f ,y/<600. Here we found 
a range 490<T c yy<600 consistent with Leggett et al., larger 
be cause of the differen t technique used. Moreover, as noted 
in IDelorme et all (2008) using the Besancon stellar population 
model dRobin et al.ll2003l) . the kinematics suggest that the target 
is a ol der member of the t hin disk (age~4 Gyr) that corresponds, 
using iBaraffe et al.l (120031) radii, to a radius R»0.9 R j up and a 
T c yy«570 K. As regards its gravity and metallicity, we noted in 
Sec. 4 that the model predictions are not consistent, so we do not 
go further in the analysis of these physical properties. 



6.3. SDSS 0207 



In IVrba et all d2004 they find a 7^=34.85+9.87 mas and 
fX tot = 156.3+1 1 .4 mas/yr for SDSS 0207, both consistent with the 
values found here. The model predictions in Figs.[2]and[3]indi- 
cates log[g]=5.0 and low or solar metallicity while the BSH06 
clear tracks indicates a metal-rich nature in J-K. 

There are t wo previous estim a tes of effective temperature 
for this object: iGolimowski et al.l (120041) . using the measured 
spectrum and estimating the bolometric corre ction using M' and 
L' photometry, find L fao/ /L = 1.5 1 ±0.87xlO' 5 : IVrba et al.ld2004 
using the K-band magnitude and the bolometric correction de- 
rived by Golimowski et al. and find L/, o //L =1.74±l.O4xlO~ 5 . 
Both these values are consistent with our result. 



6.4. ULAS 0901 

ULAS 0901 is much cooler than what expected for its type, see 
Fig. [6] In Fig. [4]the dwarf is fainter than the other T7.5s of the 
sample, in particular in My and Myy. This maybe an indication 
that U LAS 0901 is half a s ubclass later, i.e. a T8. 

In lLodieu et"all d2007l) ULAS 0901 is indicated as a high 
gravity - solar metallicity object, based on the low K-band flux. 
The model prediction of the gravity and the metallicity of this 
object in Figs.[2]and[3]are log[g]=4.5 and [Fe/H]=+0.2 (in con- 



6.5. ULAS 1018 

The model predictions in Figs. [2] and [3] are not consistent for 
ULAS 1018. In the J-K space it appears as a solar-metallicity 
low-gravity (log[g]=4.5) object, while the J-H prediction is low- 
gravity and low-metallicity ([Fe/H]=-0.5). The uncertainty in 
color and magnitude prevent us to draw a firm conclusion on 
the properties of this object. We mention here that ULAS 1018 
is indicated as a metal-poor object by Murray et al. (in p rep.) 
based on its very blue H-K, while in lLodieu et al.l was 
indicated as a possible metal-rich object, based on the high K- 
band flux and narrow Y-band flux peak. 

The bolometric luminosity and thus the temperature range 
obtained are consistent with an object of spectral type T5. 



6.6. ULAS 1150 

This object is a peculiar T6.5, with T7-T6.5 indices in the J 
and K band (H 2 Q-J = 0.087, CH 4 -J = 0.302, CH 4 -K = 0.032, 
IPinfield et al.l l2008). but T3 in the H (H 2 0-H = 0.455). Pinfield 
et al. suggest that it could be a low-gravity high-metallicity ob- 
ject, based on the K-band enhancement and the Y-band suppres- 
sion. In Fig. [2] it assumes an anomalous position, lying on the 
low-gravity high-metallicity side of the theoretical curves while 
in Figf3]the model predicts a low/solar metallicity. The large er- 
rors prevent us from making any comments on its gravity. 

The bolometric luminosity and thus the temperature range 
obtained are consistent with an object of spectral type T6.5. 



6.7. ULAS 1315 

In Fig.|6]it appears significantly cooler th an the other T7.5s. In 
the discovery paper, IPinfield et al.l d2008h . no particular indica- 
tion of peculiarity were detected. Looking at Figs. [2] and [3] it is 
easy to see that ULAS 1315 is particularly faint in Mr (it is the 
bluest in J-K), and the K-band suppression indicates high grav- 
ity, as predicted by the models. However the indication of high 
metallicity in the same plots is in contradiction with this inter- 
pretation as a metal-rich object would show an enhancement in 
the K-band flux. It is fainter than the other T7.5s also in My and 
M//, this results in an extremely cool temperature, almost 200 K 
less than what expected. This can indicate that ULAS 1315 is a 
T8 (the spectral indices H2O-J, C H4-J and CH4-H co rresponds 
to a T8, while H 2 0-H to a T7, seelEnfiddetaLlllQOll). 

6.8. ULAS 1335 

The re are two previous estima tions of the temperature of this ob- 
ject: [Burninghametal]|2008), using a solar-metallicity BTSettl 
model grid and a comparison of the H-|4 .49] color with theo - 
retical expectations, find T c// ~550-600 K: lLeggett et al.l (I2009). 
using the technique described in Sec. 6.1, find T<,yy=500-550 
K. The luminosity obtained here (L io ;=1.09+0.06xl0~ 6 L ) is 
formally consistent with both those ranges, given the large un- 
certainty in the ra d ius, b ut it is closer to the value obtained in 
iBurningham et al.l d2008l) . As in the case of the other late T 
dwarfs, ULAS 0034 and CFBDS 0059, the model predictions 
are contradictory in the different color spaces and we are there- 
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fore unable to say anything about the other parameters of this 
objects. 

7. Conclusions 

We present the parallax and proper motions of 1 1 cool T dwarfs 
taking advantage of the UKIDSS discovery image to shorten 
the time required for a precise determination. We find that the 
models do not predict the colors and absolute magnitudes of the 
coolest T dwarfs, and that, given the observed colors and ab- 
solute magnitudes, the models predict contradictory behaviours 
depending on which color is considered. 

We examine two methods for the calculation of T e ff, one of 
which derives the temperature from the luminosity, which is de- 
termined using available spectroscopic data complemented with 
model spectra scaled using measured photometry, and the other 
of which does a least-squares-like model spectrum fit to the ob- 
served spectrum. The second method is very dependent on the 
data weighting selection and so we prefer the former technique. 
The two approaches give consistent results however the former 
gives a larger range in temperature. 

The observations of these 11 objects are continuing at a 
lower frequency of 4 observations per year. At the end of the 
observing campaign all the objects in the sample will have a 
time coverage of 4 years, allowing us to obtain a more pre- 
cise and robust parallax solution. Also, enhanced photometry for 
ULAS 0901 and ULAS 0948 is expected within the year. 

The inclusion of MIR observations increases our ability to 
determine the object parameters, especially for the temperature. 
With IRAC photometry and a MIR spectrum the emitted flux 
is well constrained and the limiting factor becomes the radius, 
producing a range of ~150 K. Using a NIR spectrum and pho- 
tometry only, the uncertainty in the flux is on the order of that 
on the radius, and the temperature range obtained doubles. For 
the future we hope that the Warm-Spitzer and WISE missions 
will continue to provide MIR observations for the study of these 
objects. 
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